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We numerically study the dynamics of a discrete spring-block model introduced by Olami, Feder and Chris- 
tensen (OFC) to mimic earthquakes and investigate to which extent this simple model is able to reproduce the 
observed spatiotemporal clustering of seismicty. Following a recently proposed method to characterize such 
clustering by networks of recurrent events [Geophys. Res. Lett. 33, L1304, 2006], we find that for synthetic 
catalogs generated by the OFC model these networks have many non-trivial statistical properties. This includes 
characteristic degree distributions — very similar to what has been observed for real seismicity. There are, how- 
ever, also significant differences between the OFC model and earthquake catalogs indicating that this simple 
model is insufficient to account for certain aspects of the spatiotemporal clustering of seismicity. 
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I. INTRODUCTION 

Describing and modelling the spatiotemporal organiza- 
tion of seismicity and understanding the underlying physical 
mechanisms of earthquake triggering have been proven chal- 
lenging. Inspired by empirical regularities which include self- 
similar scaling laws like the Omori law and the Gutenberg- 
Richter law 01 01, a wealth of mechanisms and models have 
been proposed including the concepts of a critical earthquake, 
of self-organized criticality, and more generally of the seis- 
mogenic crust as a self-organized complex system requiring a 
so-called system approach (see, for example, Refs. y[ for 
a review). Yet, the origin of the observed non-trivial emer- 
gent features of earthquake occurrence is still one of the main 
unresolved problems in the field. 

Resolving this issue may require measuring the micro- 
scopic state variables — the stress and strain at every point 
within the earth along active earthquake faults — and their 
exact dynamics. This is currently impossible. However, the 
associated earthquake patterns are readily observable making 
macroscopic approaches based on the concept of spatiotem- 
poral point processes feasible, where the description of each 
earthquake is reduced to its size or magnitude, its epicenter 
and its time of occurrence. Describing the patterns of seis- 
micity may shed light on the fundamental physics since these 
patterns are emergent processes of the underlying many-body 
nonlinear system. 

Recently, such an approach has led to the identification of 
new properties of the clusterin g of seismicity in space and 
time 

SBSBIIEIIIHIIIQEIII- In particular, the 
observed spatiotemporal clustering of seismicity suggests that 
the usual mainshock/aftershock scenario — where each event 
has at most one correlated predecessor of larger magnitude 
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— is too simplistic and that the causal structure of seismic- 
ity could extend beyond immediately subsequent events. To 
quantify such correlations, a general procedure for inferring 
a plausible causal structure from clusters of localized events 
has been introduced in Refs. lfl7l [Tsh . The approach allows 
one to study the dynamical organization of spatiotemporal ac- 
tivity in terms of the topology of complex networks lfl9l I2TJI1 
and has led to the detection of unexpected statistical features 
for earthquake catalogs from California. Most importantly, 
the approach provides a new and independent estimate of the 
rupture length and its scaling with magnitude. 

In this article, we investigate to which extent the simple 
dynamics of the Olami-Feder Christensen (OFC) model 12111 
is sufficient to account for the observed spatiotemporal clus- 
tering of seismicity as characterized by the above network 
approach. The OFC model is maybe the simplest model in 
the class of self-organized critical (SOC) models which ex- 
hibit a phenomenology resembling seismicity. This includes 
the aforementioned Gutenberg-Richter law for the frequency- 
magnitude distribution I2II l22tl and the Omori law for after- 
shocks J23I l24ll as well as some statistical properties of epi- 
center locations and dynamics 1I25I1 . In addition, the OFC 
model is of special relevance in the context of SOC. Non- 
equilibrium systems are called self-organized critical if they 
reach a stationary state characterized by power laws and the 
absence of characteristic scales — without the need for fine- 
tuning an external parameter such as the temperature lEfjll27ll . 
This is typically the case for slowly driven systems with fast 
avalanche-like dissipation events. Unlike the paradigmatic 
SOC sandpile model [28], where any amount of dissipation 
is enough to introduce a scale which breaks criticality, the ex- 
istence of criticality in the OFC model with dissipation is still 
debated jHHSEllllllllll- 

We show that when the causal structure of the OFC model 
is studied with the network method mentioned above, it does 
indeed reveal some similarities to seismicity, such as the de- 
gree distributions of the network, and some aspects of the re- 
currence time and distance statistics. However, there are sev- 
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eral important differences — including the absence of a rup- 
ture length — that severely hinder its adoption as a plausible 
(sole) description of the fundamental mechanisms responsi- 
ble for the spatiotemporal correlations of seismicity. Many 
of those discrepancies seem to be closely related to the fact 
that the model is defined on a regular, discrete lattice, and 
to the existence of quasi-periodic attractors j25ll . These fea- 
tures, which are totally absent in seismicity, manifest them- 
selves quite clearly in several aspects of the recurrence net- 
work. Also, the recurrence time statistics show a characteris- 
tic time due to the fact that the system is uniformly driven. 

The outline of our paper is as follows. In Section[II] we re- 
view the OFC model, and SectionllTIlsumm arizes the network 
method to search for signs of causal structure in spatiotempo- 
ral data. Section[lV]presents the results obtained for the OFC 
model, which are compared to real seismicity. We conclude 
in SectionM 



H. OLAMI-FEDER-CHRISTENSEN MODEL 

The OFC model l2ll] is inspired by the Burridge-Knopoff 
spring -block model 13511 . and is defined on a square lattice of 
size L 2 . To each site a "tension" is assigned, initially 
chosen at random from the interval [0, z c [. The entire system 
is driven slowly, with every Zij increasing uniformly. When- 
ever a site reaches the threshold tension Zij = z c , a relaxation 
event — called "avalanche" or "earthquake" — starts. The 
originating site of the avalanche is called the epicenter. 
The dynamics of such an event is as follows: A fraction a of 
the tension at the epicenter is transferred to each of its four 
neighbors Zj±i j±i = z.;±i j±i + oiZij, and the tension at the 
site itself is reset Zij — 0. If the tension at any of the neigh- 
boring sites reaches the threshold, Zi±i,j±l > z c, the same 
toppling rules are applied again. This dynamics continues un- 
til there are no more sites in the system with Zij > z c . Then, 
the tension increases uniformly again until the next avalanche 
occurs. Without loss of generality, we set z c = 1. The total 
number of topplings during an avalanche is called the "size" of 
the avalanche s. In contrast, the total number of sites that top- 
pled (ignoring multiple topplings of the same site) is called the 
"area" of the avalanche a. The time in the system can be mea- 
sured either discretely, by the number of events ("iterations"), 
or continuously, by the cumulative tension injected in the sys- 
tem when it is driven (the so called "natural" time scale). The 
parameter a defines the level of local conservation of the sys- 
tem. For a = 0.25 the system is locally conservative and for 
a < 0.25 it is dissipative. We consider here only the case of 
open boundary conditions, i.e., the sites at the border of the 
lattice transfer tension out of the system (or to an imaginary 
neighbor), so the system is always globally nonconservative. 

After a transient regime, which is increasingly longer for 
smaller values of a, the system reaches a stationary state 
which has a distribution of event sizes resembling a power- 
law in the tail for the largest, numerically accessible system 
sizes (see Ref. ll34ll for an extensive review and the most re- 
cent results). For the values of a e {0.18; 0.2; 0.22} and 
L = 1000 we consider here, the power law tail has an unique 



exponent around —1.8 l30ll . Yet, the occurrence of events is 
not uniformly distributed within the lattice, and tend to hap- 
pen closer to the boundaries jzM l3lll . This "border effect" 
prevents the event size statistics from obeying finite-size scal- 
ing (FSS) with the size of the lattice. Only when events close 
to the border are ignored and an internal subset of the system 
is considered, FSS is realized rf3lll . Interestingly, this border 
effect does not seem to scale with system size l25ll . To ac- 
commodate for these observations, we generally restrict our 
analysis to synthetic catalogs generated by the OFC model 
which contain only those events involving sites which are all 
at a distance of at least 100 sites away from the boundary. The 
effective systems size is thus Lbuik = 800. 
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FIG. 1: (Color online) The average area (a r ) of the smallest rectan- 
gle on the lattice containing a given avalanche of spatial extent a for 
different values of the dissipation parameter a. Only avalanches con- 
tained within the bulk, Lbuik = 800, are considered (see text for de- 
tails). The solid line corresponds to a linear increase. Inset: Rescaled 
version indicating that the area scales linearly with a independent of 
a for 300 < a < 30000 (area enclosed by straight lines). 

An important aspect in seismicity is the effect of the detec- 
tion threshold: Earthquake catalogs are only complete above 
a certain magnitude value which depends on the network of 
seismometers deployed in a given area. Thus, it is crucial 
to understand how the statistical properties of seismicity vary 
with the magnitude threshold. To address this issue in the 
OFC model, we consider catalogs for different lower thresh- 
old values of event sizes s t h- Reasonable choices of these val- 
ues are limited by the particular dynamics of the OFC model 
and lattice effects. For example, it is well known that events 
of size one have their own separate dynamics and obey a dif- 
ferent statistics l34ll . This is easy to justify, since events of 
size one will occur even if there is no self-organization, i.e., 
a = 0. The effects of this particular dynamics and of the 
discrete lattice can be clearly seen in the spatial shapes of 
the avalanches as summarized in Fig. [T] The ratio between 
the average area (a r ) of the smallest rectangle which con- 
tains an avalanche of area a and a varies significantly with 
a for all considered values of a. Here the average is taken 
over all events of area a. Only for events above a certain size 
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TABLE I: Number of events N in the catalogs generated by the OFC 
model for different values of a and lower size threshold s t h- 





a — o ig 


a = 0.20 


a = 0.22 


300 


1000000 


1000000 


1000000 


500 


641933 


676324 


685237 


1000 


327124 


394235 


400339 


3000 


107550 


148311 


156189 


5000 


60041 


87102 


94968 


10000 


22205 


40208 


42991 


30000 


2547 


8432 


8909 



ath ~ 10 2 , does a r scale linearly with a (see inset of Fig. [T]) — 
independent of a. Yet, this linear scaling is limited to values 
of a < 3 x 10 4 due to the finite system size. Since area and 
size of an avalanche are basically identical for the parame- 
ters of a considered here l40ll . we only consider catalogs with 
lower size threshold s& G [300; 30000] in the following. 

To be more specific, each synthetic catalog generated by 
the OFC model contains N consecutive events or avalanches 
occurring after the statistically stationary regime has been 
reached |4lll for L = 1000 and L bu ik = 800. Each event k is 
characterized by its size Sk, its epicenter location (i,j)k and 
its time of occurrence tk measured in terms of the continuous 
natural time scale (the discrete avalanche number time scale 
of the OFC model l36ll is basically equivalent). The various 
catalogs for different values of a and s± are given in Table U 

m. NETWORK OF RECURRENCES 

We analyze the catalogs generated by the OFC model ac- 
cording to a recently proposed method which has been proven 
helpful in characterizing the spatiotemporal clustering of seis- 
micity and in detecting causal signatures between events 
lfl7[[l8ll . The essential idea is to extend the notion of a recur- 
rence to spatiotemporal point processes: An event is defined 
to be a recurrence of any previous event if it occurred closer 
in space than all intervening events. Recurrences are therefore 
record breaking events with respect to distance. By linking 
each event to its recurrences, a directed network [4211 of re- 
currences is obtained: Each event a^, with k = 1, N, is a 
vertex in the network and a directed edge from to a m ex- 
ists for k < m if a m is a recurrence of ajt, i.e., if the distance 
between (i,j) m and (i,j)k is smaller than the distance be- 
tween (i,j)k and (i,j)k' for all events as/ with k < k' < m. 
This definition assumes that the events are ordered according 
to their occurrence in time tk- Obviously, each recurrence or 
edge can be characterized by the time interval t = t m — tk 
between the two connected events k and m and by the spa- 
tial distance I between the two. Note that the mapping of the 
dynamics to a network is entirely data-driven and does not 
impose any arbitrary space and time scales other than those 
associated with the given event or earthquake catalog — in 
contrast, for example, to methods typically used to define af- 
tershocks. Comparing the statistical properties of the network 



of recurrences for a given catalog to the properties of a net- 
work obtained for a random point process without any causal 
relation between events highlights relevant parts of the under- 
lying causal dynamical process(es) generating the pattern. For 
the OFC model, such a random process or null model can be 
obtained by shuffling the entire catalog, containing all events: 
Shuffle the sizes and the epicenter locations separately while 
keeping the times of occurrence. Later the s± threshold can 
be raised, and the corresponding shuffled versions of each cat- 
alog can be obtained. For a random spatiotemporal point pro- 
cess in continuous space and time (CST), several statistical 
properties of the network of recurrences are even known ana- 
lytically Hi- 

IV. RESULTS 

In the following, we do not only compare the network prop- 
erties of the OFC catalogs with those of their random or shuf- 
fled counterparts but also with the network properties found 
for earthquake catalogs from southern California lfl7l [l8ll . 
The particular focus is on the network topology as summa- 
rized, for example, in Refs. Il9l l20ll and on distributions as- 
sociated with the temporal and spatial distances between an 
event and its recurrences. 



A. Topological properties 

1. Network growth 

An important aspect of the analysis of networks and their 
dynamics is the change in their topology with time. In the 
case of networks of recurrences, this corresponds to the sit- 
uation when a given catalog is extended to cover a longer 
time period and, thus, contains more events. In general, some 
properties of the network of recurrences depend on its size N. 
For random CST events, one such property is the average de- 
gree (k) or equivalently the average number of recurrences per 
event. It grows with the size of the network as (k) m ln(N) 
for N ^> 1 lfl8ll . While this is roughly what we find for the 
shuffled OFC catalogs, the original OFC catalogs show quite 
a different growth as can be seen in Fig. [2] Independently of 
the value of a, the average degree "saturates" quickly around 
a value close to 5, and grows very little even after the catalog 
is a couple of orders of magnitude larger. When the catalog 
is shuffled, a portion of the growth can be well represented by 
ln(N) but only up to ~ 10 4 events, after which it grows more 
slowly. 

These observations illustrate one of the main features of the 
OFC model, namely the repetitive occurrence of events origi- 
nating at the same location. Obviously, there is a finite prob- 
ability that the epicenters of two events are identical since the 
OFC model is defined on a discrete and finite lattice. If this 
happens, the cascade of recurrences of the first event ends or 
is "closed" and the number of recurrences of this event does 
not increase further (see Fig. [3]) ■ If this were exclusively due 
to the finiteness of the discrete lattice, one would expect this 
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a = 0.18 

a = 0.20 

a = 0.22 

a = 0.18, shuffled 

a = 0.20, shuffled 

a = 0.22, shuffled 




activity to return to the same epicenter and to close a given 
cascade of recurrences. Thus, (k) is larger and its growth is 
very close to the predicted behavior for random CST events 
until the effects of the broad distribution of the epicenters fre- 
quencies become important after ~ 10 4 events. 

The network growth described above differs significantly 
from what is obtained for seismicity, since real epicenters 
happen on a continuous space and therefore never occur ex- 
actly in the same place (save, of course, for limitations in pre- 
cisely locating the epicenters). Hence, the network growth 
for seismicity — at least for the given catalog sizes — does 
not exhibit the "saturation" of the average degree obtained 
for the OFC model, and grows instead continuously as (k) m 
0.841n(iV) ill. 



FIG. 2: (Color online) Mean degree (k) vs. catalog size N for s t i, = 
300 and different values of alpha, for the regular and shuffled OFC 
catalogs. N = 10 corresponds to the full catalogs — see TableU 
The solid line is of the form 0.95 In (TV) + const. 



to happen according to a geometric distribution with average 
L^ ulk . In particular, this effect would not be relevant in the 
thermodynamic limit Lbulk — ► oo. However, our simulations 
show that the repeated occurrence of the same epicenter hap- 
pens much more frequently in the OFC model. In particular, 
the OFC catalogs for s th = 300 in Table U have 79%, 87% 
and 92% closed cascades for a = 0.22,0.20,0.18, respec- 
tively i43ll . This is due to marginal synchronization of neigh- 
boring sites ll37ll and the existence of quasi-periodic patterns 
in the dynamics of the OFC model [25] causing the repeated 
triggering of epicenters and similar events. Consequently, the 
average out-degree of networks generated by the OFC model 
grows rather slowly. Moreover, the growth in (fc) basically 
saturates after ~ 10 3 events as Fig. [2] shows. This is at least 
two orders of magnitude smaller than the rough estimate of 
-^buik = 6.4 x 10 5 based on the lattice size. 
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FIG. 3: (Color online) Schematic representation of two cascades of 
recurrences of size 4. In both cases dn < dz < di < di. The 
cascade on the left is "open", i.e. it is possible for a fifth recurrence 
to occur, since di > 0. This is not the case for the cascade on the 
right. Since di = 0, this cascade is "closed" and no distance will be 
able to break the record. 



2. Degree distributions 

A more detailed aspect of the network topology is the dis- 
tribution of in- and out-degrees, P(k m ) and P(fc ollt ). In the 
context of recurrence networks, the out-degree of a vertex cor- 
responds to its number of recurrences while the in-degree is 
the number of events of which it is a recurrence. As shown in 
Fig. [4] the out-degree distribution deviates significantly from 
the random CST case, where a Poisson distribution is ex- 
pected both for the in- and out-degree distributions lfl8ll . The 
out-degree distribution seems instead to decay exponentially, 
with inclination slightly dependent on a. The in-degree dis- 
tribution is closer to a Poissonian. The overall shape of the 
distributions also does not depend on the choice of the lower 
size threshold for 300 < s th < 30000 (not shown). When the 
catalogs are shuffled, the obtained out-degree distribution is 
very close to a Poissonian (Fig. [5} as expected for the random 
CST case. 
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When the catalog is shuffled it does not entirely remove the 
effect. This is because the average frequency of a given epi- 
center does not change, and the distribution of these frequen- 
cies is quite broad, resembling a power law l25ll . Yet, the shuf- 
fling destroys the relatively short periods of quasi-periodic dy- 
namics and increases in most cases the time it takes for the 



FIG. 4: (Color online) In-degree and out-degree distributions for net- 
works generated by the OFC model for s ln = 300 and different a's 
(see TableUfor details). The solid (black) line is a Poisson distribu- 
tion with the same (k) as for a — 0.20. 

The in- and out-degree distributions obtained for the origi- 
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FIG. 5: (Color online) In-degree and out-degree distributions for net- 
works generated from shuffled OFC catalogs with s t h = 300 and dif- 
ferent a's. The solid (black) line is a Poisson distribution with the 
same (k) as for a — 0.20. 



nal OFC catalogs are basically identical to those for seismic- 
ity flUl . The only significant differences are the values of 
P(l) for the out-degree. For seismicity, it is much more likely 
to have an event with out-degree one or equivalently an event 
with a single recurrence. This is related to the finite extension 
of an earthquake, as characterized by its "rupture length", and 
the corresponding microscopic dynamics. Thus, our findings 
suggest the absence of such a length in the OFC model. This 
is indeed confirmed by other, more direct means discussed in 
Sec lIVCl Nevertheless, the close resemblance of the degree 
distributions for the OFC catalogs and earthquake catalogs 
suggests that the OFC model is able to reproduce all other 
non-trivial gross features of the spatiotemporal clustering of 
seismicity which are captured by the degree distributions. 



3. Degree-degree correlations 

The network topology canbe further characterized by the 
degree-degree correlations If2(ill . These correlations indicate 
how likely vertices of a given in- or out-degree are connected 
to vertices with another given k m or fc out . We also consider the 
correlations between km and fc out of the same vertex, the sin- 
gle vertex degree-degree correlations. To capture both types 
of correlations, we measure the average k m or fc out (either of 
the neighbors or of the same vertex) where the average is taken 
over all vertices with a given km and fc out . For the OFC cat- 
alogs, Figs. [6] and [7j show the absence of pronounced non- 
trivial (single vertex) degree-degree correlations. There are 
only small qualitative differences between the shuffled and un- 
shuffled catalogs [44] indicating that the degree-degree corre- 
lations are not particularly sensitive to the given form of spa- 
tiotemporal clustering. This is in agreement with findings for 
seismicity 1181 . Interestingly, the single vertex degree-degree 
correlation function (fc out ) (km) increases with fc; n (Fig. [7} — 



both for the shuffled and the original OFC catalogs. This 
is different from shuffled earthquake catalogs and especially 
from random CST processes which show a slight decrease re- 
lated to the age of a vertex lll8ll . The deviating behavior of the 
(shuffled) OFC catalogs is due to the closing of cascades of 
recurrences mentioned above which largely reduces the corre- 
lation between the age of a vertex and its in- and out-degree. 
Instead, vertices with large k m correspond to sites with low 
activity which do not belong to the quasi-periodic attractors. 
These sites do not experience the aforementioned closing of 
the cascade of recurrences, and thus have a larger out-degree. 
The shuffling of the catalog does not alter this since it does 
not affect the closing of cascades — only the time it takes to 
close. This reasoning is confirmed by the right panel of Fig. [7] 
which shows the single vertex degree correlation, when only 
open cascades of recurrences are considered, and the increase 
is less evident — and disappears completely for the shuffled 
case. 
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FIG. 6: (Color online) Degree-degree correlations for a = 0.2 and 
Sth = 300. The curves on the bottom (green) are for the original 
catalog, and on top (blue) for the shuffled catalog. Note that the offset 
between the curves for the original OFC catalog and the shuffled 
version is simply due to different mean degrees. 
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FIG. 7: (Color online) Single vertex degree-degree correlations for 
a — 0.2 and so, = 300. The left panel shows the results for the 
original catalogs, and the right panel shows the results for the open 
cascade of recurrences only. The bottom curves (green) are for the 
unmodified catalogs, and the top ones (blue) are for the shuffled cat- 
alog. 
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TABLE II: Clustering coefficients (Eq.Q}) of the recurrence net- 
works generated by the OFC model for different values of a and 
lower size threshold stb = 300. The values of C excluding events 
with closed cascade of recurrences and for the shuffled OFC catalogs 
are shown for comparison. 



a 


C 


C (open) 


C (shuffled) 


0.18 


0.4153(8) 


0.2064(8) 


0.15411(8) 


0.20 


0.3910(3) 


0.1970(6) 


0.15891(8) 


0.22 


0.3470(3) 


0.1891(4) 


0.14650(7) 



4. Clustering coefficient 

The clustering coefficient l20ll is another quantity that char- 
acterizes the network topology. It quantifies how likely neigh- 
bors of the same vertex are connected with each other. In the 
present context, this refers to the probability that recurrences 
of the same event are also recurrences of each other. 

For all vertices i with out-degree larger than one, the local 
clustering coefficient C, is given by the ratio of existing links 
Ei between its fc° ut recurrences to a maximum possible num- 
ber of such links, ^k° ut (k° at — 1). The clustering coefficient 
C of the network is defined as the average over all vertices i 
with out-degree larger than one 

c={Ci)i = \kr(lr-i)) i ' (1) 

The values of the clustering coefficient for the different 
OFC catalogs are given in table [U] The original catalogs 
generate networks which are significantly more clustered than 
those generated by the shuffled catalogs. This effect decreases 
slightly for larger values of a. A substantial amount of clus- 
tering seems to be related to the quasi-periodic attractors of 
the dynamics generating vertices with closed cascades of re- 
currences. This follows directly from the observation that the 
clustering is much less pronounced for those vertices with 
open cascades of recurrences (Table iHli. However, the latter 
events still retain a higher clustering than the networks for the 
shuffled OFC catalogs. 

While a direct quantitative comparison with seismicity is 
not possible due to different catalog sizes and different mean 
degrees, the qualitative features are the same. In both cases, 
the catalogs show a much higher clustering compared to their 
respective shuffled counterpart. 



B. Temporal distances of recurrences 

The probability density function (PDF) p s ' h (t) for the time 
intervals or waiting times t associated with the recurrences 
or the edges of the network is another important characteris- 
tic. Fig. [8] shows this PDF for different OFC catalogs with 
a = 0.2. Despite some variation with s±, all PDFs decay 
approximately as 1/t up to the largest possible waiting time, 



determined by the finite time span of the catalog, with a super- 
posed peak at a characteristic time t c . As the inset of Fig. [8] 
shows, this characteristic time separates two slightly differ- 
ent regimes. While the PDFs decay almost exactly as 1/t 
for t > t c , the decay is rather 1/t ' 9 for t < t c . Note also 
that t c does not depend on s±. Such a characteristic time is 
expected since the activity shows quasi-periodic behavior as 
discussed in Section llV A II Moreover, t c depends on the dis- 
sipation parameter a, as Fig. [9] shows. More specifically (as 
the inset of Fig. [9] shows), the characteristic time seems to 
be exactly 1 — 4a, which is the natural period of the system: 
For periodic boundary conditions and relatively small a, the 
system exhibits only trivial, perfectly synchronized topplings 
of size one and the activity returns to each site with period 
1 — 4a |38|. It is even possible to identify secondary charac- 
teristic recurrence times in the inset of Fig. |9j corresponding 
to integer multiples of 1 — 4a, again emphasizing the quasi- 
periodic behavior Ir45tl . 
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FIG. 8: (Color online) Distribution of the waiting times between 
events and their recurrences for different threshold sizes s t h, for 
a = 0.2. The solid line on the left is a power law with exponent 
—0.9 and on the right one with exponent — 1. Inset: Rescaled distri- 
butions with iV taken from table U 

For the shuffled OFC catalogs, p Slh (t) is shown in Fig. [TOl 
As expected, the peak at the characteristic time scale t c is 
barely noticeable since the quasi-periodic behavior is basi- 
cally absent in the shuffled catalog. Yet, all other qualitative 
features are preserved. In particular, the PDF is very different 
from the behavior of p Slh (t) in the random CST case where a 
flat regime for small t and 1/t decay for large t exist lfl8ll . In 
the latter case, the transition time t* between the two regimes 
scales with the inverse of the rate of events t* ~ T/N. We 
do not observe such a scaling for the shuffled OFC catalogs as 
follows from the inset of Fig. [10] While the data do collapse 
for large arguments, there is no good collapse for small argu- 
ments. This is due to the fact that there is a non-negligible 
probability that the time interval between subsequent events 
in the OFC model is arbitrarily small i36ll39ll — in sharp con- 
trast to the assumption of a Poisson process in the random 
CST case. 
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FIG. 9: (Color online) Distribution of the waiting times between 
events and their recurrences for Sth = 300 and varying a. Inset: 
Rescaled distribution, zoomed around the characteristic time region. 
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FIG. 10: (Color online) Distribution of waiting times between events 
and their recurrences for shuffled OFC catalogs with different thresh- 
old sizes Sth and a = 0.2. The solid line on the left is a power law 
with exponent —0.9 and on the right one with exponent — 1. Inset: 
Rescaled distribution with TV taken from tableU 



For seismicity t* — and indeed p Slh (i) — does not depend 
on N (corresponding to a varying magnitude threshold) and 
there is also no equivalent to t c IU8I1 . These observations 
clearly prove that the OFC model has certain features which 
are incompatible with seismicity and which strongly influence 
the generated spatio temporal clustering of events. 



Here, it is assumed that the recurrences for a given event are 
ordered according to their time of occurrence such that the z-th 
recurrence is followed by the i + 1-th recurrence for all ranks 
i. The corresponding PDFs are shown in Fig.Qj] For the OFC 
catalogs, panel (a) indicates that the PDFs vary with i. While 
they closely resemble a power law with exponent —0.85 in 
an intermediate regime for small i, deviations occur for larger 
i's. This is independent of a and s t h (not shown). Independent 
of i, there is also a peak around ti/ti+i = 1. As Fig. [TTI(c) 
shows, there are only few discernible differences between the 
original and the shuffled catalogs: The peak is slightly more 
pronounced in the original catalog implying that subsequent 
recurrences are more likely to be separated by a very short 
time interval only, and the occurrence of smaller ratios for 
larger values of i is more likely for the shuffled catalogs. A 
third and more pronounced but expected difference is that the 
kink related to tc around ti/ti+i ~ 10~ 10 in the original cat- 
alog is absent in the shuffled catalog. 

Even these small differences in the PDFs disappear if only 
open cascade of recurrences are considered, as shown in 
Fig.[TT](b) and (d). In all cases the PDFs depend very weakly 
on the value of a, with all mentioned properties equally 
present in the range of a's studied. 




io 8 




- i = 1 
» i = 2 
» i = 3 
■ i = 5 


io 6 




- i = 7 


io 4 






io 2 

10° 


. (b) 






io" 12 io 8 io" 4 


10° 


io 8 






io 6 






io 4 






io 2 

10° 


(d) "''■3a* 





10 



10 10 

ti/ti+1 



10 



FIG. 11: Distribution of the waiting time ratios for subsequent recur- 
rences, for a = 0.2 and s t h = 300, for (a) the regular catalog, (b) 
only the "open" recurrences, (c) the shuffled catalog and (d) only the 
"open" recurrences for the shuffled catalog. In all plots the solid line 
is a power-law with inclination —0.85. 



1. Temporal hierarchy of subsequent recurrences 

Another important aspect of the network of recurrences is 
the possible existence of hierarchies of recurrences. This can 
be captured by the ratios of the waiting times ti/U+i for sub- 
sequent recurrences belonging to the same originating event. 



For seismicity, the ratios U/ti+i for subsequent recurrences 
revealed a different non-trivial characteristic Their dis- 
tribution showed two power-law regimes — in contrast to a 
single power law for the shuffled catalog. The spatiotempo- 
ral correlations responsible for these results do not seem to be 
present in the OFC model, as the above analysis shows. 
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C. Spatial distances of recurrences 

Exactly like the waiting times of recurrences, the spatial 
distances between an event and its recurrences can be obtained 
for the OFC model for further characterization and compar- 
ison with seismicity. The PDF of these distances p s,h (l) is 
shown in Fig. Q~2] for a = 0.2 and different size thresholds. 
The results for other values of a are qualitatively similar. The 
overall shape of the distributions is rather broad with cut-offs 
at short and long distances due to the discrete and finite lat- 
tice, respectively, and varies with s±. Yet, for arguments less 
than the long distance cut-off two different regimes seem to 
be generally present. The inset of Fig.Q~2]shows an attempted 
scaling collapse according to the ansatz 



p s HI) 



(l/si)-PF(l/s 



6 ) 

th) 



(2) 



with 5 = 0.3 and /3 = 1.5. If successful for arguments sig- 
nificantly less than the respective finite size cut-off, it would 
imply the existence of a characteristic distance dependent on 
s,h with I* ~ Sfa 3 — qualitatively similar to seismicity and the 
random CST case as discussed below. However, as the inset of 
Fig.Q~2]clearly shows, the data do not collapse as intended and 
a collapse does not seem to be possible at all. When only open 
cascades of recurrences are considered, as shown in Fig. Qj] 
the PDFs are slightly different but there is no convincing col- 
lapse either. It does, however, show again that the statistics 
of those events that have closed cascades and predominantly 
belong to quasi-periodic attractors, are significantly different 
from the "open" ones. 
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FIG. 12: (Color online) Distribution of recurrence distances for a = 
0.2 and different threshold values s t h- Inset: Distributions rescaled 
according to Eq.{2} with S — 0.3 and f3 = 1.5. Recurrences of 
distance I — are not shown, and were not used in the normalization. 



For the shuffled OFC catalog the corresponding results for 
p Slh (l) are shown in Fig. [14] In this case, a collapse is in- 
deed possible. Here, the collapse was done using the number 
of events TV of the respective catalog given in Table U Due 
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FIG. 13: (Color online) Distribution of recurrence distances for a = 
0.2 and different threshold values s t h, but only considering "open" 
cascades of recurrences. Inset: Distributions rescaled according to 
Eq.{2} with 5 = 0.3 and (3 = 1.35. 



to the absence of non-trivial correlations in the shuffled cata- 

-0.8 



-1.8 



logs, this is equivalent to the scaling with s t h since TV oc s, 
as follows from the distribution of event sizes P(s) ~ 
These findings can be directly compared to the random CST 
case for which I* oc TV -7 , where 7 = 1/D and D is the 
underlying spatial dimension flal . However, for the shuffled 
OFC catalog we find I* oc TV -018 which would imply D sw 5. 
This is in sharp contrast to the fact that the OFC model is 
defined on a 2D lattice. This disagreement can be attributed 
to the assumption of continuous space made to derive the re- 
sults in Ref. fl8tl . On a discrete lattice, when the character- 
istic distance approaches the fundamental lattice distance, the 
identification of the exponent with 1 /D cannot be expected. 
Despite this mismatch, the random CST process correctly pre- 
dicts the power-law decay with exponent 1 for large distances 
of the shuffled OFC catalogs. Note that this decay clearly dis- 
tinguishes the shuffled OFC catalogs from the original OFC 
catalogs. 

The equivalence of the scaling of p Slh (l) with s± and with 
TV for the shuffled OFC catalogs discussed above is confirmed 
by Fig. [15] (c). There, the variation of p Slh (l) with the time 
length At of the catalog is explicitly shown for s t h = 300. 
Note that the x and y axis are rescaled with At 1 , where At is 
the time length indicated in the legend. Not only the shuffled 
catalog varies with time (and hence TV) for small distances 
but also the original catalog — even if only open cascades 
of recurrences are considered (see Fig. [15] (a) and (b)). In 
the latter case, the different curves can be collapsed using the 
same ansatz as for the shuffled catalog with 7 = 0.37, while 
the time dependence of the unmodified catalog seems to be 
more elaborate (it is shown unsealed in Fig. [T5l ). Moreover, 
comparing Figs. [15] (b) and the inset of [13] allows to test the 
scaling relation TV 7 oc s^ 0,87 in the case of open cascades 
of recurrences. Indeed, we have 0.8 x 0.37 ~ 0.3 assuming 
TV oc At. This would imply that the dependence of p s ' h (l) 
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to identify I* with the rupture length of an earthquake — its 
spatial extent — even though it is explicitly neglected in the 
point process description. In particular, it directly relates the 
underlying microscopic dynamics to the statistical properties 
of the network of recurrences — one of the main achieve- 
ments for seismicity. As shown above, the OFC model does 
not successfully reproduce these features. While the charac- 
teristic distance does depend on the lower threshold size, no 
simple scaling ansatz can be identified (Eq.(f2|), even when 
only "open" cascades of recurrences are considered. More 
importantly, the characteristic distance decreases with time in 
both cases. 



1. Spatial hierarchy of subsequent recurrences 



FIG. 14: (Color online) Distribution of recurrence distances for the 
shuffled OFC catalog with a = 0.2 and different threshold values 
s t h- Inset: Distributions rescaled according to the number of events 
(see Table 13 as indicated in the legend. 



on s t h is simply due to the variation in N. Yet, the fact that 
p s,h (I) does not obey scaling suggests that the situation is more 
complicated and that there is a non-trivial dependence on s±. 
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FIG. 15: (Color online) Distribution of recurrence distances for a = 
0.2 and s t h = 300, for an increasingly larger potion of the catalog 
At, re-scaled as indicated in the legend, for (a) the regular catalog 
(7 = 0, i.e. no scaling), (b) only the "open" recurrences (7 = 0.37) 
and (c) the shuffled catalog (7 = 0.18). 



The spatial hierarchy of recurrences can be captured by 
the distribution of distance ratios U+i/U for subsequent re- 
currences belonging to the same originating event. Here, it 
is assumed that the recurrences for a given event are ordered 
according to their time of occurrence such that the i-th recur- 
rence is followed by the i + 1-th recurrence. The correspond- 
ing PDFs are shown in Fig. [16] for a = 0.2 and s± = 300 
but the PDFs seem to be independent of the specific values 
of a and s t h in the range we studied. We have also included 
the case i = for which we use Iq = \/2Lb u ik- Note that l\ 
corresponds to the distances between subsequent events and 
that in this case there is a cut-off at large arguments due to the 
finite lattice size. For all z > the PDFs shown in Fig. [16] 
seem to be independent of i. For the original OFC catalog 
(Fig.fToTa)), the PDFs increase for small arguments and then 
become flat for larger ones, before increasing again as the ra- 
tio approaches 1, The same behavior is observed when only 
open cascades of recurrences are considered (Fig.[T6jb)). 

The situation is different for the shuffled OFC catalogs. As 
shown in Fig.fTSTc), all PDFs increase as a power-law with ex- 
ponent 1, This is similar to the random CST model which pre- 
dicts p(x) = Dx D ~ x , independent of i, with 2: = h+i/k IU8I1 . 
For the shuffled OFC catalogs, this would imply a spatial di- 
mension D = 2 — exactly as expected. The same behavior 
is observed when only open cascades of recurrences are con- 
sidered (Fig. [T6i d)). The comparison of the OFC catalogs and 
their shuffled counterparts suggests that the distribution of ra- 
tios between subsequent recurrence distances indeed captures 
non-trivial aspects of the spatiotemporal dynamics. 

All cases are however very different from seismicity, where 
the ratios all decay as a power-law with exponent ^0.6. 



V. DISCUSSION & CONCLUSION 



For seismicity, p Slh (l) shows a characteristic distance I*, af- 
ter which the distribution decays as a power-law with an ex- 
ponent close to one fl^l . The value of I* depends only on the 
lower magnitude threshold and, in particular, does not change 
with time (or N) for fixed m. More specifically the scaling of 
this characteristic distance obeys the ansatz l*(m) oc io - 45m , 
where m is the lower threshold magnitude. This allows one 



We have studied the network of recurrences of the OFC 
model, with the objective of characterizing its causal struc- 
ture and comparing it to seismicity. We have found that, in 
agreement with previous results I23ll24ll25]l . the model does 
exhibit some non-trivial features observed in seismicity, but 
it comes short of a complete qualitative description of several 
robust features of spatiotemporal clustering captured by the 
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FIG. 16: Distribution of ratios between subsequent recurrence dis- 
tances for a — 0.2 and s t h = 300, for (a) the regular catalog, (b) 
only the "open" recurrences, (c) the shuffled catalog and (d) only the 
"open" recurrences for the shuffled catalog. The solid line is a power 
law with exponent 1. 



feature, as both the time and distance distributions vary with 
magnitude threshold and time, respectively. The second fea- 
ture is also not reproduced. While there is a dependence of 
the distance distribution on the size of the events, there is no 
scaling which would allow to identify a rupture length. One 
third, and more subtle, characteristic of the recurrences which 
is not reproduced by the OFC model is the hierarchy of re- 
currences, given by the distribution of ratios of time and dis- 
tance between consecutive recurrences, which in the case of 
seismicity deviates significantly from the random case. More- 
over, the whole analysis exposes another crucial aspect of the 
OFC model, which is the elevated occurrence of epicenters 
at distance zero from one another, tightly connected to quasi- 
periodicity and well distinct from seismicity. 

Obviously, the OFC model was not intended to be a com- 
plete description of seismicity. It was mainly proposed as 
a conceptual model and a possible origin of the Gutenberg- 
Richter law from a SOC point of view. Also, like the 
Burridge-Knopoff spring-block model, it aims to model the 
dynamics of a single fault, and not the complex interplay of 
different fault structures as observed in reality. Despite these 
limitations, the OFC model is capable of showing a rich dy- 
namics, reproducing even some more subtle features of the 
spatiotemporal clustering of earthquakes. 



used method. More specifically, the model reproduces cen- 
tral topological features of the recurrences network, which are 
the degree distributions and high clustering, as well as the ab- 
sence of non-trivial degree correlations and some other par- 
tial aspects of the recurrence distance and time distributions. 
However, the two most important properties of the dynamics 
of seismicity identified by this recurrence approach are, firstly, 
that the time (distance) distributions are independent of mag- 
nitude thresholds (time), and secondly, that the distance statis- 
tics provides an independent assessment of the rupture length 
of earthquakes. The OFC model fails to reproduce the first 
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